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Abstract. - Twisted rapid passage is a type of non-adiabatic rapid passage that gives rise to 
controllable quantum interference effects that were first observed experimentally in 2003. We show 
that twisted rapid passage sweeps can be used to implement a universal set of quantum gates that 
operate with high-fidelity. For each gate in the universal set, sweep parameter values are provided 
which simulations indicate will yield a quantum gate with error probability P e < 10~ 4 . Note that 
all gates in this universal set are driven by a single type of control field (twisted rapid passage), 
and the error probability for each gate falls below the rough-and-ready estimate for the accuracy 
threshold P a ~ 10~ 4 . The simulations suggest that the universal gate set produced by twisted 
rapid passage shows promise for use in a fault-tolerant scheme for quantum computing. 
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Introduction. — The physical context for our discus- 
sion is the accuracy threshold theorem [1-8] which estab- 
lished that a quantum computation of arbitrary duration 
could be done, with arbitrarily small error probability, 
in the presence of noise, and using imperfect quantum 
gates, under the following conditions. (1) The computa- 
tional data is protected by a sufficiently layered concate- 
nated quantum error correcting code. (2) Fault-tolerant 
protocols for quantum computation, error correction, and 
measurement are used. (3) A universal set of unencoded 
quantum gates is available with the property that each 
gate in the set has an error probability P e that falls be- 
low a value P a known as the accuracy threshold. The 
value of the threshold is model-dependent, though for 
many, P a ~ 1CF 4 has become a rough-and-ready estimate. 
Thus gates are anticipated to be approaching the accura- 
cies needed for fault-tolerant quantum computing when 
P e < 10~ 4 . One of the principal challenges facing the 
field of quantum computing is finding a way to implement 
a universal set of unencoded quantum gates for which all 
gate error probabilities satisfy P e < 10~ 4 . 

In this Letter numerical simulation results are presented 
which suggest that a class of non-adiabatic rapid pas- 
sage sweeps, first realized experimentally in 1991 [9], and 
known as twisted rapid passage (TRP), should be capa- 



ble of implementing a universal set of unencoded quantum 
gates Q u that operate non-adiabatically, and with gate er- 
ror probabilities satisfying P e < 10~ 4 . Q u consists of the 
one-qubit Hadamard and NOT gates, together with vari- 
ants of the one-qubit 7r/8 and phase gates, and the two- 
qubit controllcd-phase gate. The universality of Q u was 
demonstrated in Ref. [10]. This level of gate accuracy is 
due to controllable quantum interference effects that arise 
during a TRP sweep [11,12], and which were observed us- 
ing NMR in 2003 [13]. To find sweep parameter values 
that yield such high-performance gates, it proved neces- 
sary: (i) to combine the simulations with an optimiza- 
tion procedure that searches for minima of P e [10,12]; and 
(ii) for the modified controllcd-phase gate, to also apply 
the symmetrized evolution of Ref. [14] . 

The outline of this Letter is as follows. We begin with 
a summary of the essential properties of TRP. This is fol- 
lowed by a discussion of how the simulation and optimiza- 
tion are done, and how symmetrized evolution is incorpo- 
rated into the two-qubit dynamics. We then present our 
simulation results for the different gates in Q u . We close 
with a discussion of our results and of future work. 

Preliminaries. — To introduce TRP [11,12], we con- 
sider a single qubit interacting with an external control 
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field F(t) via the Zeeman interaction H z (t) = cr ■ F(t), 
where the o~i are the Pauli matrices (i = x, y, z). TRP is a 
generalization of adiabatic rapid passage (ARP). In ARP, 
the field F(t) is slowly inverted over a time To such that 
F(i) = atz + 6x. In TRP, the control field is allowed to 
twist in the x-y plane with time-varying azimuthal angle 
4>(i), while simultaneously undergoing inversion along the 
z-axis: 

F(i) = atz + b cos (f>(t)x + b sin <j>(t)y. (1) 

Here -T /2 < t < T /2, and the TRP inversion can be 
non-adiabatic. 

Controllable Quantum Interference. As shown in 
Rcf. [12], the qubit undergoes resonance when at — 
?>4(t)/2 = 0. For polynomial twist with <j> n (t) = (2/n)Bt n , 
this condition has n — 1 roots, though only real-valued 
roots correspond to resonance. Ref. [11] showed that for 
n > 3, the qubit undergoes resonance multiple times dur- 
ing a single TRP sweep: (i) for all n > 3, when B > 0; and 
(ii) for odd n > 3, when B < 0. For the remainder of this 
Letter we restrict ourselves to B > 0, and to quartic twist 
for which n — 4 in 4> n (t). For quartic twist, the qubit 
passes through resonance at the times t = 0,±^a/nB 
[11]. Thus the time separating the qubit resonances can 
be controlled through variation of the sweep parameters 
B and a. Ref. [11] showed that these multiple resonances 
have a strong influence on the qubit transition probabil- 
ity, allowing transitions to be strongly enhanced or sup- 
pressed through a small variation of the sweep parameters. 
Ref. [15] calculated the qubit transition amplitude to all 
orders in the non-adiabatic coupling. The result found 
there can be re-expressed as the following diagrammatic 
series: 



r_(t) = 




Lower (upper) lines correspond to propagation in the neg- 
ative (positive) energy level, and the vertical lines corre- 
spond to transitions between the two energy levels. The 
calculation sums the probability amplitudes for all inter- 
fering alternatives that allow the qubit to end up in the 
positive energy level at time t given that it was initially 
in the negative energy level. As we have seen, varying the 
TRP sweep parameters varies the time separting the res- 
onances. This in turn changes the value of each diagram 
in eq. (2), and thus alters the interference between alter- 
natives in the quantum superposition, It is the sensitivity 
of the individual alternatives/diagrams to the time sepa- 
ration of the resonances that allows TRP to manipulate 
this quantum interference. Zwanziger et al. [13] observed 
these interference effects in the transition probability us- 
ing NMR and found quantitative agreement between the- 
ory and experiment. It is the link between the TRP sweep 
parameters and this quantum interference that we believe 
makes it possible for TRP to drive highly accurate non- 
adiabatic one- and two-qubit gates. 



Simulation and Optimization. A detailed presentation 
of our simulation and optimization protocols appears in 
Refs. [10, 12]. We can only give a brief sketch of that 
presentation here. As is well-known, the Schrodinger dy- 
namics applies a unitary transformation U to an initial 
quantum state \tp) which is driven by the system Hamil- 
tonian H(t). The Hamiltonian (see below) is assumed to 
contain terms that Zceman-couple each qubit to the TRP 
control field F(t). Assigning values to the TRP sweep 
parameters (a,b,B,To) determines H(t), which then deter- 
mines the actual unitary transformation U a applied. The 
task is to find sweep parameter values that produce a U a 
that approximates a target gate Ut sufficiently closely that 
its error probability (defined below) satisfies P e < 10~ 4 . 
In the following, the target gate Ut will be one of the 
gates in the universal set Q u . Since Q u contains only one- 
and two-qubit gates, our simulations will only involve one- 
and two-qubit systems. For the one-qubit simulations, the 
Hamiltonian H\(t) is the Zccman Hamiltonian introduced 
earlier. Rcf. [10] showed that it can be written in the 
following dimcnsionless form: 

H\{t) = (1/A) {-tg z - cos 04 (t)^ - sin 04 (t)^} . (3) 

Here: r = (a/b)t; X — ha/b 2 ; and for quartic twist, 
4>4 (t) = (?74/2A)t 4 with r]4 = %Bb 2 /a? . For the two-qubit 
simulations, the Hamiltonian H 2 (t) contains terms that 
Zeeman-couple each qubit to the TRP control field, and 
an Ising interaction that couples the two qubits. Note 
that alternative two-qubit interactions can easily be con- 
sidered, though we focus on the Ising interaction here. 
The energy-levels for the resulting Hamiltonian contain a 
resonance-frequency degeneracy that was found to spoil 
gate performance. Specifically, the resonance-frequency 
for transitions between the ground- and first-excited states 
(Ei <-> E%) is the same as that for transitions between 
the second- and third-excited states (£3 «-> E4). To re- 
move this degeneracy a term c^E^t)) (E±(t)\ was added 
to H2(t). Combining all these remarks, one arrives at the 
following dimcnsionless two-qubit Hamiltonian [10]: 

«2(r) = 

[-(di + d 2 )/2 + r/X]al - (d 3 /A)[cos <£ 4 ct£ + singer*] 
+ [-d 2 /2 + T/X]a 2 z - (1/A) [cos (j) A a 2 x + sin 4 a 2 ] 
-(nd4/2)al<j 2 z + C4 \E 4 (t))(E 4 (t)\. (4) 

Here: (1) h = h^iB rf /2, uji = jiB Q , and i = 1, 2; (2) r = 
(a/b2)t, A = ha/b 2 ,, and 774 = hBb 2 /a 3 ; and (3) di = (u>\ — 
w 2 )6 2 /a, d 2 = (A/a)6 2 , d 3 = 61/62, and d A = (J/a)b 2 , 
where A is a detuning parameter [10]. 

The numerical simulation assigns values to the TRP 
sweep parameters and then integrates the Schrodinger 
equation to obtain the unitary transformation U a pro- 
duced by the sweep. To assess how closely U a approxi- 
mates the target gate U t , it proves useful to introduce the 
positive operator P = (U^ — U})(U a — U t ). Given U a , 
Ut, and an initial state \ip), one can work out the error 
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probability P e (t/j) for the TRP final state \tp a ) = U a \ip), 
relative to the target final state \ipt) = Ut\ip). The gate 
error probability P e is defined to be the worst-case value 
of P e (tp): P e = max| V)) P e (tp). Rcf. [12] showed that P e 
satisfies the bound P e < Tr P, where the RHS is the trace 
of the positivie operator P introduced above. Once U a is 
known, Tr P is easily evaluated, and so it makes a conve- 
nient proxy for P e , which is harder to calculate. To find 
TRP sweep parameter values that yield highly accurate 
non-adiabatic quantum gates, it proved necessary to com- 
bine the numerical simulations with function minimiza- 
tion algorithms that search for sweep parameter values 
that minimize the Tr P upper bound [16]. The multi- 
dimensional downhill simplex method was used for the 
one-qubit gates, while simulated annealing was used for 
the two-qubit modified controllcd-phase gate. This pro- 
duced the one-qubit gate results that will be presented 
below. However, for the modified controllcd-phase gate, 
simulated annealing was only able to find parameter val- 
ues that gave P e < 1.27 x 10~ 3 [10]. To further improve 
the performance of this two-qubit gate, it proved necessary 
to incorporate the symmetrized evolution of Ref. [14] to 
obtain a modified controlled-phase gate with P e < 10~ 4 . 
We now briefly describe how symmetrized evolution is in- 
corporated into our simulations. 

Symmetrized Evolution and TRP. Ref. [14] intro- 
duced a unitary group-symmetrization procedure that 
yields an effective dynamics that is invariant under the 
action of a finite group Q. We incorporate this group- 
symmetrization into a TRP sweep by identifying the group 
Q with a finite symmetry group of the target gate Ut, 
and then applying the procedure of Ref. [14] to filter out 
the ^-noninvariant part of the TRP dynamics. As the Q- 
noninvariant dynamics is manifestly bad dynamics relative 
to Ut, group-symmetrized TRP yields a better approxima- 
tion to Ut- We briefly describe the group-symmetrization 
procedure, and then show how it can be incorporated into 
a TRP sweep. 

Consider a quantum system Q with time-independent 
Hamiltonian H and Hilbcrt space TC. The problem is to 
provide Q with an effective dynamics that is invariant un- 
der a finite group G, even when H itself is not ^-invariant. 
This symmetrized dynamics manifests as a Q -invariant ef- 
fective propagator U that evolves the system state over a 
time t. Let {pi = p{gi)} be a unitary representation of 
Q on H, and let \G\ denote the order of Q. The proce- 
dure begins by partitioning the time-interval (0,t) into N 
subintervals of duration A/jjv = t/N , and then further par- 
titioning each subinterval into \Q\ smaller intervals of du- 
ration Stpj = A£jv/|Cf|. Let SUn = cxp [— (i/f^St^H] de- 
note the _ff-generated propagator for a time- interval Stpj, 
and assume that the time to apply each pi £ Q is neg- 
ligible compared to 8tN (bang-bang limit [17]). In each 
subinterval, the following sequence of transformations is 
applied: U{At N ) = ]Tl=i p\&U N pi. Rcf. [14] showed 
that: (i) {/(Atpj) — > exp[— (i / h) At n H] as TV — > oo, where 



H = (1/IGI) Ei='i PiHpi; (ii) H is ^-invariant {[H, Pi ) = 
for all pi £ Q); and (iii) the propagator U over (0,i) is 
U = cxp[— (i/h)tH], which is C/-invariant due to the Q- 
invariance of H . The end result is an effective propagator 
U that is (J-invariant as desired. 

This procedure can be generalized to allow for a time- 
dependent Hamiltonian H(t). To do this, the time in- 
terval (0,t) must be divided into sufficiently small subin- 
tervals that H(t) is effectively constant in each. Within 
each subinterval, the above time-independent argument 
is applied, yielding a ^-symmetrized propagator for that 
subinterval. Combining the effective propagators for each 
of the subintervals then gives the full propagator U = 
T{cxp(—i/h J Q cItH(t))], where T indicates a time-ordered 

exponential, and H(t) = p\H(t)p t . 

For our two-qubit simulations, the target gate is 
the modified controlled-phase gate V cp = (1/2) [(J 1 + 
crl)I 2 — (I 1 — o-\)o-1] which is invariant under the group 
Q = {J 1 / 2 , a*, o%, a\o-l}. Thus \Q\ = 4, and we set 
pi = I 1 / 2 ,...,/^ = o-^cr 2 . Switching over to di- 
mensionless time, we partition the sweep time-interval 
(—To/2, To/2) into sufficiently small subintervals that 
our two-qubit Hamiltonian HzIt) is effectively constant 
within each. We then apply the time-independent sym- 
mctrization procedure to each subinterval with the V cp 
symmetry group acting as Q. Combining the effec- 
tive propagators for each of the subintervals as above 
gives the ^/-symmetrized propagator for the full TRP 
sweep U = T[exp(-i/h) f^°J o 2 /2 dTH(T))}, with H(t) = 

(^/^)J2i=iPi^( T )Pi- We shall see that ^-symmetrized 
TRP yields an approximation to V cp with P e < 10 -4 . 

Simulation Results. 

One-qubit gates. Operator expressions for the tar- 
get gates are: (i) Hadamard — Uh = (l/V2)(o~z + o- x )', 
(ii) NOT — U not = a x ; (iii) modified n/8 — V„/ s = 
cos(7r/8)cr a; — sin(7r/8)f7j,; and (iv) modified Phase — V p = 
(l/v / 2)(c a: — o~ y ). The gate fidelity is calculated using 
T n = (l/2 n )Re[Tr(W a U t )] 1 where n is the number of 
qubits acted on by the gate. This fidelity is especially 
convenient as it is related to our Tr P upper bound: 
T n = 1 - (l/2 n+1 )TrP [10]. Finally, the connection be- 
tween the TRP experimental and theoretical parameters 
is given in Refs. [10], [11], and [12] for superconducting, 
NMR, and atom-based qubits, respectively 

A study of the TRP-implementation of these one-qubit 
gates was first reported in Rcf. [12]. It proves useful 
to rcparamcterize the TRP sweep parameters A — > A* 
and 774 — ► 77I, where A = A*cxp[-(A* - A°)/A°] and 
?/4 = 77I cxp [—(77! —wD/wl]- For each one-qubit gate in Q u , 
the fixed-point of the rcparametcrization (A , rfl) is given 
by the optimum sweep parameter values found in Rcf. [12]. 
Table 1 presents the values for the optimum sweep param- 
eters (A*, 77!) that produced our best results for Tr P for 
each of the gates in Q u . The fixed-point (A , for each 
gate is also listed. In all one-qubit simulations, the dimen- 
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Table 1: Simulation results for the one-qubit gates in Q u . The gate error probability P e satisfies P e < TrP. 



Gate 


A* 


vt 




TrP 


A° 






u h 


5.85 


2.93 x 10 


-4 


9.30 x 10" 


(j 


5.8511 


2.9280 x 10 


-4 


Unot 


7.32 


2.93 x 10 


-4 


1.12 x 10" 


5 


7.3205 


2.9277 x 10 


-4 


U-k/8 


6.02 


8.15 x 10 


-4 


3.55 x 10" 


5 


6.0150 


8.1464 x 10 


-4 


u p 


5.98 


3.81 x 10 


-4 


8.70 x 10- 


5 


5.9750 


3.8060 x 10 


-4 



Table 2: Variation of TrP for the Hadamard gate when the TRP sweep parameters are altered slightly from their optimum 
values. Variation of Tr P for the other one-qubit gates in Q u is comparable to that of the Hadamard gate and so corresponding 
Tables for these other gates are not shown. 



A* 


vl 




TrP 




??4 


A* 


TrP 




5.85 


2.92 x 10 


-4 


2.24 x 10" 


5 


2.93 x 10" 4 


5.84 


1.24 x 10" 


-5 




2.93 x 10 


-4 


9.30 x 10" 


-6 




5.85 


9.30 x 10" 


6 




2.94 x 10 


-4 


6.06 x 10" 


•5 




5.86 


1.12 x 10" 


5 



sionless inversion time was r = 80.0. Since P e < Tr P, we 
see that P e < 10~ 4 for all one-qubit gates in Q u . The Tr P 
values yield the following gate fidelities: (i) Hadamard — 
T h = 0.9999 98; (ii) NOT — T not = 0.9999 97; (iii) mod- 
ified tt/8—^/s = 0.9999 91; and (iv) modified Phase— 
J- p = 0.9999 78. In Table 2 we show how Tr P varies when 
A* and rjl are varied slightly from their optimum values for 
the Hadamard gate. Similar behavior occurs for the other 
one-qubit gates in Q u , and so in the interests of brevity, 
corresponding Tables for these other gates are not shown. 
We see that gate performance is a slowly- varying function 
of the sweep parameters (A* , 774 ) . 

Modified controlled-phase gate. We complete the uni- 
versal set Q u by presenting our simulation results for 
the ^-symmetrized TRP implementation of the modified 
controlled-phase gate V cp . In the two-qubit computational 
basis (eigenstates of cr^crf), V cp = diag(l, 1,-1,1). TRP 
implementation of V cp without symmetrized evolution was 
reported in Ref. [10]. The results presented there are su- 
perceded by the Cy-symmetrized TRP results presented be- 
low. For purposes of later discussion, note that the pa- 
rameters appearing in T-L^t) fall into two sets. The first 
consists of the TRP sweep parameters (A, 774, to), while 
the second set (a, di,. . . , d&) consists of parameters for 
degeneracy-breaking, detuning, and coupling. We parti- 
tioned the TRP sweep into N seq = 2500 pulse sequences, 
with each sequence based on the four element symmetry 
group for V cp introduced earlier. 

For the modified controlled-phase gate V cp , gate perfor- 
mance was not found to be very sensitive to small vari- 
ations of the TRP sweep parameters. Instead, for V cp 
without symmetrized evolution [10], gate performance was 
most sensitive to C4, d\, and d±. However, when sym- 
metrized evolution was incorporated into the TRP sweep, 
d\ ceased to be a critical parameter. Thus it only proved 
necessary to reparameterize C4 — ► c| and d± — ► d|, where 
c 4 = c|exp[-(c4-c^)/c^] and rf 4 = ci|exp[-(d| -d%)/d%\. 
Simulations incorporating symmetrized evolution deter- 



mined the reparameterization fixed-point to be c° = 2.173, 
e?4 = 0.8347. The (optimized) parameter values A = 5.04, 
r] 4 = 3.0x 10~ 4 , r = 120, di = 99.3, d 2 = 0.0, d 3 = -0.41, 
d\ = 0.835, and c\ = 2.17 produced a gate U a for which 
TrP = 8.87 x 10" 5 , gate fidelity T cv = 0.9999 89, and 
gate error probability satisfying P e < 8.87 x 10~ 5 . We see 
that by adding symmetrized evolution to a TRP sweep 
we obtain an approximation to V cp with P e < 10 -4 . In 
Table 3 wc show how Tr P varies when the parameters 
C4 and d^ vary slightly from their optimum values. Sen- 
sitivity of gate performance to the remaining parameters 
is comparable to that of c% and d\ and so corresponding 
tables are not shown. We see that gate performance is a 
slowly- varying function of the parameters c\ and d%, as 
well as of (A, 774, To) and (di, d?, d 3 ). 

Discussion. — We have presented simulation results 
which suggest that TRP sweeps should be capable of im- 
plementing a universal set of quantum gates Q u that oper- 
ate non-adiabatically and with gate error probability satis- 
fying P e < 10~ 4 . Wc note that all gates in the universal set 
Q u are driven by a single type of control field (TRP), and 
that the gate error probability for all gates in Q u falls be- 
low the rough-and-ready estimate for the accuracy thresh- 
old P a ~ 10~ 4 . The simulation results presented in this 
Letter suggest that the universal quantum gate set Q u pro- 
duced by TRP shows promise for use in a fault-tolerant 
scheme for quantum computing. Refs. [10-12] have shown 
how TRP sweeps can be applied to NMR, atomic, and 
superconducting qubits. It should also be possible to ap- 
ply them to spin-based qubits in quantum dots using a 
magnetic field since the same Zeeman-coupling acts as 
with NMR qubits. Although we have studied a number of 
forms of polynomial twist, as well as periodic twist [18], 
we have found that quartic twist provides best all-around 
performance when it comes to making the gates in Q u - 
At present we do not have arguments that explain why 
quartic twist works better than other forms of twist. We 
are currently working to develop a theory of the optimum 
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Table 3: Variation of TrP for the modified controlled-phase gate when c\ and d% are altered slightly from their optimum values. 



r* 


d% 


TrP 


dl 


cl 


TrP 




2.17 


0.834 


8.77 x 10- & 


0.835 


2.16 


8.15 x 10" 


-5 




0.835 


8.77 x 10- 5 




2.17 


8.77 x 10" 


-5 




0.836 


8.77 x 10~ 5 




2.18 


8.44 x 10" 


-5 



twist profile in an effort to understand this question. 
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